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ABSTRACT 

Relationships  between  the  real  and  imaginary  parts  of  an  impedance  are 
studied.  Given  the  real  and  imaginary  parts  in  a  finite  frequency  band,  a  technique 
is  developed  for  finding  these  parts  outside  the  band.  Conditions  under  which  the 
real  part  remains  positive  for  all  u  are  given. 

Use  of  a  negative  impedance  converter  for  realization  of  a  particular  im¬ 
pedance  arising  from  impedance  loading  studies  is  discussed.  A  method  of  com¬ 
pensating  for  collector  capacitance  and  extending  operation  of  a  transistorized 
NIC  into  the  VHF  region  is  given. 
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I 

INTRODUCTION 

This  report  is  devoted  to  those  aspects  of  the  impedance  loading  problem  which 
are  concerned  with  network  theory  rather  than  with  field  theory. 

In  seeking  to  control  the  radar  cross-section  of  a  conducting  body  by  impedance 
loading,  the  scattering  behavior  of  the  body  is  first  determined  as  afurction  of 
loading  impedance,  either  lumped  (as  for  a  thin  wire)  or  distributed  (as,  for  exam¬ 
ple,  with  a  sphere).  The  equations  are  then  solved  to  give  the  impedance,  as  a 
function  of  frequency,  necessary  to  obtain  the  desired  cross-section  modification. 
The  network-theory  pr^lern  is  that  of  synthesizing  and  realizing  this  impedance 
from  numerical  data  giving  its  real  and  imaginary  parts  over  some  band  of  fre¬ 
quencies. 

Network  synthesis  actually  involves  two  separate  problems,  approximation  and 
realization.  The  given  data  must  first  be  approximated  by  a  rational  function  of 
frequency  chosen  from  an  appropriate  class  of  rational  functions  (e.g. ,  positive- 
real  functions).  A  network  realizing  this  rational  function  is  then  obtained  by  any 
of  several  fairly  standard  techniques. 

The  approximation  problem  arising  from  the  impedance-loading  study'  is  a 
rather  unusual  one.  Both  the  real  and  imaginary  parts  of  the  imj)edance  are  speci¬ 
fied  over  some  relatively  large,  but  finite,  band  of  frequencies;  the  more  common 
situation  is  for  some  one  quantity,  such  as  magnitude,  to  be  specified  over  the 
entire  frequency  range. 

The  unasurl  nature  cf  the  problem  has  prompted  a  theoretical  study  of  the 
relations  between  parts  of  impedances  specified  over  a  finite  band.  A  number  of 
interesting  and  useful  results  have  been  obtained,  and  are  reported  in  the  next 
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sections.  These  results  are  expected  to  form  the  basis  for  techniques  of  handling  a 
large  class  of  network  problems  arising  in  connection  with  reactive  loading. 

An  experimental  program,  aimed  at  more  immediate  practical  results  in  a 
specific  case,  has  also  been  undertaken.  This  program  has  attempted  to  obtain  an 
active-network  realization  of  an  impedance  whose  behavior  is  approximated  by  a 
pure  negative  reactance;  such  an  impedance  was  found  by  Chang  and  Senior  (1967)  to 
be  required  to  minimize  the  cross  section  of  a  sphere  loaded  by  a  slot  parallel  to 
the  incident  wave. 

The  principal  effort  toward  realization  of  a  negative  reactance  has  been  in 
the  area  of  negative  impedance  converters  (NIC’s).  Considerable  success  has  been 
achieved  in  compensating  for  high-frequency  effects' and  extending  the  useful  range 
of  operation  of  the  NIC  into  the  VHF  range.  These  results  are  discussed  in  Section  T 
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n 

THEORETICAL  STUDIES 

The  principal  theoretical  problem  arising  from  the  study  of  impedance  loading 
can  be  broadly  formulated  in  the  following  terms: 

Given  the  real  part,  R  (u) ,  and  the  imaginary  part,  X  (u>) ,  of  an  impedance 
on  some  frequency  interval,  what  can  be  said  concerning  the 

behavior  of  these  functions  outside  the  band  ? 

A  number  of  more  specific  questions  can  be  based  on  this  general  statement  of 
the  problem.  For  example: 

Does  there  exist  a  function  Z  (p),  analytic  in  the  right-half  p-piane,  whose 
real  and  imaginary  parts  are  equal  to  F. (u)  and  X(w),  respectively,  for  p  -  ja, 

US  <  U)  <  U)  ? 

1  —  —  2 

If  such  a  function  exists,  is  it  unique  V 

Assuming  that  Z  (p)  exists,  what  procedure  can  be  used  to  find  it,  at  least 
for  p  =  j  u  ? 

What  are  the  conditions  under  which  Z  (p)  is  positive-real:  that  is,  under  what 
conditions  will  Re  £z  (j  wj]  >  0  for  all  «  be  satisfied  ? 

All  the  questions  have  been  considered,  some  of  them  by  several  different 
methods.  All  of  them  have  been  answered  at  least  in  part;  Z{p)  does  exist,  it  is 
unique,  Z(ju)  can  be  found  fcr  all  us,  and  some  necessary  (though  not  sufficient) 
conditions  for  Z  (p)  to  be  positive-real  are  known. 

2. 1  General  Method  of  Attack 

The  most  powerful  analytic  tool  in  this  investigation  has  been  the  Hilbert  Trans¬ 
form,  which  relates  the  real  and  imaginary  parts  of  a  function  analytic  in  the  closed 
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right-half  plane.  If  Z  (p)  is  such  a  function,  then 

|  z(p, 


_  p  - 
C  o 


dp  =  0 


(2.1) 


where  C  is  the  contour  shown  in  Fig.  2-1.  Letting  the  radius  of  the  small  semi¬ 
circle  approach  zero  and  that  of  the  large  semicircle  approach  infinity,  and  sepa¬ 
rating  real  and  imaginary  parts,  we  obtain 


life) 


X{«) 


-if 

•if 


r 

X(X)dX 

J  -co 

X-u) 

a  oo 

R(X)dX 

L 

X  -u 

*  R(co) 


(2.2a) 


(2.2b) 


These  integrals  are  improper,  and  must  be  evaluated  in  the  Cauchy  principal-value 
sense,  that  is. 


p  CD 

*-co 


dX  = 


lim 


[r:  •  *  c  “] 


(2.3) 


The  basic  Hilbert  Transform  pair,  (2. 2a)  and  (2. 2b),  maybe  manipulated  into 
a  number  of  other  useful  forms.  In  particular,  if  the  integrals  are  expressed  as 
the  sura  of  two  Integrals,  one  from  -co  to  0  ami  the  other  from  0  to  ©,  one 
obtains,  after  a  suitable  change  c£  variable 


Rfo)  =  -  |  f 


XX  (X)  dX 

2  2 
X  -w 


R(oo) 


<2. 4a) 
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A  oo 
2  id  f 

*  J0 


R(X)  dX 

.2  2 

X  -id 


(2.4b) 


where  *  se  has  been  made  of  the  fact  that  R  (id)  is  an  even  function  and  X  (w)  an 
odd  function. 

The  Hilb '  rt  Transform  need  not  be  applied  directly  to  Z  (p) .  An  especially 
useful  set  of  relations  is  obtained  by  applying  the  transform  to  the  function 
Z(p)  //*  +  p2  .  Of  the  four  equations  thus  obtained,  the  two  of  greatest  value  in 
this  work  are 


■  1  [-f 1  ■ 

1  LJ  o  (x2-u2)/i-: 


X((d)  = 


R  (X)  dX 
(X2-«Ml  -XT 


X  R(X)  d  X 

(x2-id2)  yj 

X(X)  dX 

(x2-u2)  ix 


=-h| ,  o <w< 

-ij 

(2. 

,  0  <  id 


(2.5b) 


A  simplification  which  has  been  used  in  much  of  this  investigation  is  that  of 
assuming  R(id)  and  X(id)  known  in  the  band  0  <  id  <  1  rather  than  for  u^uKidg. 
This  assumption  results  in  considerable  simplification  of  the  mathematical  manip¬ 
ulations  without  any  significant  loss  of  generality,  since  the  results  obtained  for  the 
"  low  -pass”  case  can  be  readily  extended  to  the  "band-pass”  case. 

The  fundamental  problem  of  this  investigation  has  been  the  solution  of  the 
integral  equations,  (2.  5a)  and  (2. 5b),  for  the  unknown  parts  of  R(td)  and  X(id). 
Although  closed  forms  of  solution  have  not  been  found,  a  number  of  numerical 
techniques  have  been  developed;  these  are  discussed  in  succeeding  sections  of  this 
report. 
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2.2  Bounds  on  Function  Behavior 

It  is  reasonable  to  ask  what  restrictions  are  placed  on  R  (u)  and  X  (u)  in 
(0, 1)  under  the  assumption,  necessary  for  physical  realizability  of  a  passive  net¬ 
work,  that  R  (to)  be  non-negative  for  all  w  .  From  (2. 4b) 


X(u)  = 


_  2u  r1 

'  io 


R  (X)  d  X 

.2  2 
X  -u 


P00  R(X)  dX 

Ji  x2-.2 


(2.6) 


With  R(u)  >  0,  the  second  integral  in  (2.6)  is  a  non -negative  ,  non-decreasing 
function  of  u  for  0  <  u  <  1 .  Therefore 

X(u)  >  —  f  K^!  i*  ,  0  <  u  <  1  (2.7) 

'’Jo  X2-u 

>  A.  ["it!  f  R*X)dXl,  0  <  u  <  1  .  (2.8) 

dU  -  da  j0  X2-U2J 

By  adding  and  subtracting  R  (u)  inside  the  integral  we  get 

— —  f!  Sf#  =  hi  f1  .  I  SHln(B  ,2.9) 

'  Jo  x2-.2  1  Jo  X2-.2  1  W 

where  the  first  integral  on  the  right  is  now  proper.  This  result  indicates  that  the 
bounds  of  (2.  7)  and  (2. 8)  are  rather  weak,  since  the  function  on  the  right  has  a 
logarithmic  singularity,  approaching  -  co  as  w  — >  1  -  .  For  example,  for 

R(u)  =  — —-2  ,  X(u>)  =  — — «  ,  0  <  u  <  1  (2.10) 

1+u  1+u 
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we  have 


1  +  u 


2  - 


>  - 


u) 


1+u 


1  A  1  .  + 

2  +  75  ,n  \rrsy 


(2.11) 


and  also  that  the  slope  of  the  left-hand  term  of  (2. 11)  must  exceed  that  of  the  right- 
hand  term.  These  two  functions  are  plotted  in  Fig.  2-2,  which  shows  the  weakness 
of  the  inequality  near  w  =  1 . 

It  should  be  pointed  out  that  inequalities  (2.  7)  and  (2. 8)  become  much  stronger 
if  R  (1)  =  0,  but  of  course  this  is  a  special  case. 

A  stronger  inequality  can  be  obtained  by  observing  that  the  second  integral  on 
the  right  side  of  (2. 5a)  is  non-negative  for  0  <  to  <  1 .  Thus  we  have 


RH  >  -  f1  . 

J0  (X-u>fl-X 


0  <  u  <  1 


(2.12) 


or  more  generally 


R  (id)  >  -  4  ^o2-u/‘  |  ° 


XX  (X)dX 


r2  2  J~2~ '  2’ 

(a  -u  nfu  -X 


0  ^  It}  u) 


(2.13) 


Using  the  same  pair  of  functions  as  the  previous  example; 
1 


R(w)  = 


X(u))  =  - 


1  +  u 


u 


1  +  u 


(2.14) 


the  inequality  (2. 13)  is  shown  in  Fig.  2-3  for  =  1 ,  and  -  2 .  It 

will  be  noted  that  as  u  increases,  the  lower  bound  approaches  the  function  itself, 

o 

as  one  would  expect.  Letting  u- — >00  in  (2.13)  and  comparing  with  (2.4a),  we 
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see  that  the  lower  bound  actually  approaches  R  (u)  -  R(oo);  this  makes  no  dif¬ 
ference  in  the  example  or  indeed  in  most  cases  of  practical  interest. 

If  the  right  side  of  (2. 13)  goes  negative  anywhere  on  (0,  (indicating  that 
R(co)  is  non-zero),  a  better  lower  bound  can  be  obtained  by  adding  a  suitable  con¬ 
stant  to  make  the-  minimum  value  of  the  lower  bound  zero. 

2. 3  Fourier  Methods 

One  way  of  approaching  the  problem  of  finding  R(y)  and  X(u)  outside  the  band 
in  which  they  are  given  is  by  means  of  the  change  of  variable 

u  =  (2.15) 


which  maps  the  right-half  p-planc  into  the  unit  disc  in  the  y-  plane.  If  Z  (p)  is 
analytic  for  Re(p)  >  0,  then  Z(u)  will  be  analytic  for  ]u)|  <  l.  Thus,  a  power 
series  expansion  around  the  origin 


(2.16) 


converges  for  1  u>)  <  1 .  On  die  unit  circle,  let 

„  -  ,-50 


Then 


Z(€‘j*) 


+  A.  c"^  +  -r 


=  Ao  A  S  c  "  *  V 


=  (Ao  +  AjCOS0  +  A0  cos  2(5  +  ...  ) 


-  j (A^  sin0  +  A0  sin 20  +  ...) 


(2.17) 


(2.18) 
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From  (2.15)  and  (2.17)  we  have 

0 

u  =  tan  i?  ,  -jr  <  0  <  tt 

30  that 

R~{0)  =  R(tan  ~)  =  A  +  A.cos0  +  A  cos 20  +  ... 

Z  O  1  z 

X*(0)  =  X(tan|-)  =  B1sin0+  B2sin20  +  ... 

=  -A1sin0  -  A2sin20  +  ... 


(2.19) 


(2.20) 


(2.21) 


We  thus  have  the  Wiener- Lee  criterion  for  physical  realizability 

A  +  B  -  0  ,  n  =  1,  2,  ...  .  (2.22) 

n  n 

If  R(w)  and  X(cj)  are  given  for  ue  (0, 1),  then  R"(0)  and  X  "(0)  are  known 

for  -  <  0  <  — ■  .  (The  extension  to  negative  values  of  0  comes  from  the 

z  z 

fact  that  R  and  X  are  respectively  even  and  odd. )  One  possible  method  of  solutioi 
is  to  find  a  Fourier  Series  which  gives  the  best  fit  (in  the  least-square  sense)  to 

i»r  «*» 

R  (0)  and  X(0),  subject  to  the  constraint  (2. 22).  This  method  was  used  by 
Redheffer  (1947);  unfortunately,  the  equations  are  rather  difficult  to  generate,  and 
the  determination  of  R"  and  X  '"  for  |0  J  >  requires  the  summation  of  a 
Fourier  Series. 

The  present  investigation  made  use  of  a  different  approach.  Assume  that  the 
range  <  0  <  jr  is  divided  into  a  equal  intervals,  and  that  R*(0)  and  X  :;:(0) 
in  each  o£  these  intervals  may  be  approximated  by  a  constant.  1  bus  we  have 
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R  '”'{0)  =  Rk 

X*(0)  =  x k 


•  ■ 


ri  =  £  +  k  — 
*k  2  2n 


The  Fourier  coefficients  are  then  gi  /en  by 


(2.24) 


=  fjj  ^R^jOcosmpd 0  +  Z7  J** 

*  'IC  X<'(*,sinm',d<,  J 


R=i:(0)cosm0d0  +  f  ^ 

k  =  i  <!0k_. 


k  =  1  *  0. 
r 


Rj_cosm0  d0 


X,  sin  m 0  d0 


(2.25) 


(2.26) 


Applying  the  Wiener  -  Lee  criterion  we  obtain,  after  evaluating  the  right-hand 
integrals  in  (2. 25)  and  (2. 26) 


(0)  cosm0  +  X^{0)  sinmjy  d0 


k-l 


sin m 0  -sinm0  cosm0,  -cosm0, 

k  k  - 1  k  r) 


(2.27) 


For  m  =  1, 2,  ....  2t,  (2. 27)  yields  a  set  of  2n  simultaneous  linear  equations  for 

the  2n  unknowns  R  ,  — ,  R  ,  X  ,  — ,  X  .  The  coefficient  matrix  is  independent 
l  n  1  n 

of  RT  and  X "  ,  although  it  depends  on  n  . 

Aside  from  the  approximate  nature  of  the  solution,  this  method  has  the  dis¬ 
advantage  that  R  (u)  and  X(u)  are  found  as  the  real  and  imaginary  parts  of  an 
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impedance  whose  behavior  is  restricted.  From  (2. 15)  ami  (2. 16)  we  have 


Since  all  the  poles  of  Z  (p)  are  at  p  =  -1 ,  Z  (p)  lacks  the  approximating  power  of 
a  general,  nth  order  function. 

The  computational  difficulties  of  this  method  are  avoided  fey  using  an  alternative 
Fourier  method  discussed  in  the  next  section. 


2.4  Alternative  Fourier  Method 

Equations  (2. 5a)  and  (2. 5b)  may  be  written  in  the  form 


f(ca)  = 


g(u)  = 


«  Rfcj) 

if- 7 


f 


AX  (X)  dA 


r 


AR(X)'dA 


1  (X2-.2)^2-! 


2d 


X(u)  j 

j*1  R(X)dA 

f°  X(A)  dA 

(  2*  J 

fl  -u 

'°  ,.2  2sl  v2* 

(A  -u  )|i  -A 

11  (X2-u2,1 

1  2  ‘ 
[A  - 1 

0  <  t>  <  1 

(2.2Sa) 

0  <«<  1 

(2.29b) 


where  f(u)  and  g(u)  are  known  functions  for  0  <  u  <  1 . 

Since  X(A)  is  an  odd  function,  it  may  be  written  as  A E (A),- where  E(A)  is 
even.  This  substitution  makes  the  right-hand  sides  erf  (2. 29a)  and  {2. 29b)  formally 
identical,  so  that  it  will  suffice  to  consider  only  one  erf  them. 
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Making  the  change  of  variable 
X  =  sec  6 
in  {2. 29a)  we  have 

'/2 
0 


{2.30} 


ffc>) 


•j: 


R  *  (0)  d  9 

2  2 
1  -  u  cos“  0 


(2.31) 


where 


R^Ce)  =  R  {sec  9) 


{2.32) 


Note  that  if  R^{0)  >  0  for  0  <  6  <  :;/2,  the  right-hand  side  of  (2.31),  and 
all  its  partial  derivatives  with  respect  to  is,  are  non-negative  for  0  <_  y  <  1 . 
Therefore,  if  f  (is)  is  expanded  in  a  power  series  about  the  origin 


oo 


f  (jj)  -  7  K  iii 
n 

n  =  0 


2n 


(2.33) 


the  K  *s  must  all  be  non-negative. 
n  2  2-1 

The  term  (1  -u  cos  “  9)  appearing  on  the  right  side  of  {2. 31 )  may  be  ex¬ 
panded  by  the  binomial  theorem  into  an  infinite  series  which  converges  absolutely 
for  u  <  1.  Thus  we  have 


co 

-E 


K 


^n 


n  =  0 


r/2  -  ^ 

L  R'<e)Z> 


4n 

cos  &d  0  . 


{2.34) 


“0 


Since  the  series  converges  absolutely.,  we  may  interchange  the  order  of  integration 
and  summation  and  equate  coefficients  of  like  powers,  obtaining 
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where 


R’(?)  = 


R{sec  |-) 


(2.41) 


Since  R’{0)  is  a ft  even  function  it  can  be  expanded  in  a  Fourier  cosine  series 


R’{0 


A 


cos  m0 


(2.42) 


m  =  0 


Making  use  of  the  trigonometric  identity 

2n  0  1  (  ,  2n. 

cos  2  =  -fD  It^O  5  COS'° 


-J>?  +  <2„n}| 


(2.43) 


ami  substituting  (2.42)  and  (2.43)  into  (2.4C)  we  obtain 

n- 1 


K  -  1 

n  2n+ 1 


f  T?  Acosm  2  (2n)  cos (n  -  j)  0  -  (2na)l 

J  0  “o  l  j=0  J  J 


►d0 
(2.44) 


From  the  orthogonality  relation 


|  cosm0  cos(n- j)  0  d0  =  ^ 

J  0 


-  ,  m  =n-j  =  0 
~  ,  m  =n- j / 0 


(2.45) 


V  0  ,  m  ^  n-j 


we  obtain  finally 


Kn  ‘  22n+l  n-j  ) 


Z7  *  .<’n> 


(2.46) 
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This  relation  is  solved  recursively  to  obtain 
.,2n+ 1 


^n  t 


K 


It  would  appear  that  by  making  use  of  the  relation 
n 

P 

_» 

k  =  0 


cosntf  =  y  '  b  cos2k  ^ 
4— nk  2 


where 


b  =2 
nn 


2n-  1 


.  _  /  lVn-k  2n  /n+k-l\  _2k-l  ,  , 

bnk  -  H)  (  n-k-l  )  2  •  ■ 


in  (2. 42),  one  could  obtain  an  expression  of  the  form 
oo 


R’(0)  =  T?  B,  cos2k 


k  =  0 

which  corresponds  to 


oo 


R(w)  =  B,  /  u 


2k 


k  =  0 


(2.47) 


(2.48) 


(2.49) 

(2.50) 


(2.51) 


(2.52) 


However,  such  a  form,  since  it  is  a  power-series  expansion  about  infinity,  will  be 
meaningful  only  if  Z(p)  has  no  singularities  outside  the  unit  circle.  Even  for 
cases  satisfying  this  condition,  serious  convergence  difficulties  were  encountered, 
most  probably  due  to  the  truncation  of  Fq.  (2.42). 
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This  difficulty  is  easily  circumvented  by  summing  the  Fourier  Series,  (2.42), 
directly,  and  then  changing  the  abcissa  for  0  to  w  by  use  of  the  relation 


w  =  sec 


£ 

2 


(2.53) 


Results  obtained  by  applying  this  technique  to  known  functions  have  been  ex¬ 
cellent.  The  real  and  imaginary  parts  of  a  known  impedance,  calculated  on  (1,  oo) 
using  a  seven-term  series,  agreed  with  the  actual  values  to  three  significant  digits. 

The  error  due  to  truncation  can  be  estimated  rather  easily.  We  note  from 
(2.  46)  and  (2.47)  that  the  relationship  between  ne  first  (N+l)  A^'s  and  K^’s  is 
independent  of  the  higher-order  terms.  Thus,  truncating  the  Fourier  Series  (2.42) 
at  the  N  th  term  is  equivalent  to  truncating  the  power  series  (2. 34)  at  the  N  th  term. 
A  measure  of  the  error  is  thus  given  by 


e  (u)  =  f(u) 


-£  v2"  ■  Z 

n=0  n  =  N+ 1 


K  u 
n 


2n 


(2.54) 


where  e(u)  is  positive  for  all  u,  since  the  K^'s  are  all  positive. 
To  avoid  difficulty  with  singularities,  we  define 

XX  (A)  dA 


F(u)  ^  j  R(u)  + 


f?r 

» 0 


(2. 55) 


2  J  z 
■u  )  11 -X 


and 


e:!(u)  =  F(u)  -  /V  K 

n 

n  =0 


where  e  *  (u>)  >  0  ,  0  <  u  <  x  . 


2n  I  21 

|1  -U) 


(2.56) 
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We  now  define  the  error  by 


E 


Defining 

I 


e*(w)  du> 


F{u)  dw 


F(u»)  du 


and  carrying  out  the  integration  in  the  numerator  of  (2. 57)  we  obtain 

£  yn  (2n- 1)  (2n-3)  ...  (3)  (1) 

2Zj_  n  (2n+2)(2n)  ...(4)(2) 


(2.57) 


(2.58)  j 


(2.59) 


2.5  Direct  Computer  Solution 

A  number  of  attempts  were  made  to  obtain  a  direct  solution  by  means  of  the 

computer.  Given  Rq(u)  and  Xq(w)  for  me  (w^,  ui^),  the  objective  is  to  develop 

a  computer  program  which  will  determine  either  the  coefficients  of  a  positive-real 

rational  function  or  the  element  values  of  a  passive  network,  such  that  the  real  and 

imaginary  parts  of  the  corresponding  impedance  approximate,  in  some  sense,  R  (to) 

o 

and  X  (u). 
o 

The  most  general  approach  is  by  way  of  a  positive-real  rational  function.  The 
method  which  was  used  may  be  summarized  by  the  following  steps: 
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1)  As  a  starting  point,  pick  the  coefficients  of  a  rational  function,  Z(p),  of 
degree  n: 


Z(p)  = 


a  P  +  ...  +  a  p  +  a 
n  _ 1  o 

b  p  +  . . .  +  b.  p  +  b 
n  1  o 


<2. 60) 


The  program  will  reject  Z  (p)  if  it  is  not  positive-real . 

2)  Evaluate  the  real  and  imaginary  parts  of  Z{p),  R(u.)  and  X(u.)  respec¬ 
tively,  at  the  m+ 1  frequencies  u. ,  where 


u  =o)  <  u  <  ...  <  w  =  u, 

a  o  1  mb 


3)  Compute  the  weighted  error 


(2.61) 


wi  (R(“i’  ■  v^)2  +  (x<“i>  -  vv)2 


(2. 62) 


where  the  w.'s  are  the  weights. 

4)  Determine  the  partial  derivatives  of  E  with  respect  to  the  coefficients 

ao*  ar  "**  '  an’  V  br  ‘ ‘ (2.  60)  and  compute  the  gradient  direction  in 
the  (2n  +  2)  -  dimensional  coefficient  space. 

5)  Move  in  the  negative-gradient  direction  in  coefficient  space  until  one  of  the 
following  happens: 

a)  The  error  reaches  a  minimum.  In  this  case,  compute  the  new  gradient 
and  proceed  as  above. 

b)  The  function  ceases  to  be  positive-real.  In  this  case,  retreat  along  the 
gradient  line  back  to  the  positive-real  region,  compute  the  gradient 
again,  and  proceed  as  above. 
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6)  When  a  point  is  reached  from  which  no  further  progress  can  be  made  (the 
error  cannot  be  reduced  further  while  remaining  in  the  positive-real  region)  the 
search  terminates  and  results  are  printed. 

In  practice,  this  procedure  was  found  to  be  very  slow.  This  seemed  to  be  due 
to  the  foe:;  that  the  boundaries  of  the  positive-real  region  in  the  coefficient  space  are 
quite  irregular,  requiring  frequent  backtracking  and  gradient  computation.  In 
addition,  the  process  of  checking  the  function  for  positive-real  character  is  quite 
time  consuming.  In  addition,  since  the  error  is  a  very  complicated  function  of  the 
coefficients,  it  was  felt  that  local  minima  would  make  the  process  quite  dependent 
on  the  initial  values  of  the  coefficient.  This  line  of  investigation  was  therefore  not 
pursued  farther. 

A  simpler,  although  less  general,  procedure  is  to  choose  a  network  configuration 
and  adjust  the  element  values.  The  condition  for  physical  realizability  is  now  much 
simpler;  all  the  element  values  must  be  non-negative.  The  procedure  which  was 
used  was  a  gradient  search  (steepest  descent  method)  similar  to  that  outlined  above. 
The  most  significant  difference  was  that  the  boundary  of  the  region  of  physical 
realizability  is  much  more  easily  defined  in  the  present  case.  Whenever  movement 
in  the  negative-gradient  direction  caused  one  or  more  element  values  to  become 
negative,  they  were  simply  set  back  to  zero  before  the  next  step  was  taken.  Tnus, 
the  search  technique  was  actually  a  modified  gradient  search,  since  the  direction 
of  search  deviated  from  the  gradient  line  near  the  boundaries  at  the  positive  - 
coefficient  region. 

This  “direct  synthesis”  technique  had  much  better  convergence  properties  than 
the  rational  function  approximation.  The  improvement  appeared  to  be  due  to  the 
simple  nature  of  the  constraints  in  this  case  (element  values  non -negative)  compared 
to  those  of  the  previous  case  (function  positive-real).  A  number  of  difficulties 
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remained,  however.  The  method  is  obviously  dependent  on  the  configuration  chosen 
for  the  network;  a  ladder  network  was  used  in  this  study  for  its  computational 
simplicity,  but  this  may  not  be  the  best  choice  from  a  practical  standpoint.  The 
problem  of  local  minima  and  dependence  on  initial  values  still  remained.  Further, 
it  became  apparent  that,  due  to  the  high  dimensionality  of  the  problem  and  the 
extremely  nonlinear  dependence  of  the  error  function  on  the  element  values,  the 
gradient  search  technique  was  not  a  sufficiently  powerful  method. 

A  quadratic  search  was  then  tried.  In  contrast  to  the  gradient  search,  which 
essentially  linearizes  the  function  to  be  minimized,  the  quadratic  search  retains 
another  term  of  the  power  series  expansion.  This  technique  will  minimize  a 
quadratic  form  in  a  single  step.  Unfortunately,  the  problem  at  hand  did  not  involve 
a  quadratic  form  but  a  much  more  complicated  function,  to  which  the  method  was 
totally  unsuited.  In  fact,  the  procedure  was  unstable,  making  large  jumps  in 
element-value  space  without  getting  an\r  closer  to  the  minimum  point. 

No  further  effort  was  devoted  to  this  area  of  direct  computer  solution.  It  was 
quite  obvious  that  much  more  powerful  search  techniques  would  be  required 
involving  much  longer  computation  times.  Judging  from  the  results  obtained  thus 
far,  prospects  for  success  are  poor,  even  with  the  use  of  more  powerful  search 
techniques. 
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EXPERIMENTAL  PROGRAM 

Since  the  theoretical  studies  discussed  above  are  of  rather  long-range  appli¬ 
cability,  an  experimental  study  aimed  at  more  immediate  results  was  conducted  in 
parallel.  Efforts  were  directed  toward  finding  a  network  realization  of  the  imped¬ 
ance  shown  in  Fig.  3-1.  This  impedance  was  determined  by  Chang  and  Senior  (1967) 
as  the  required  loading  impedance  for  cross-section  reduction  of  a  sphere  loaded  by 
a  circumferential  slot  in  the  plane  erf  incidence. 

In  the  region  of  interest,  near  ka  =  1,  the  desired  impedance  can  be  well 
approximated  by  a  negative  reactance,  which  suggests  realization  by  an  active  net¬ 
work.  Unfortunately,  there  exists  no  really  coherent  theory  for  the  realization  of 
active  driving-point  impedances;  the  task  was  therefore  one  of  choosing  from  a 
number  of  available  alternatives  those  which  seemed  to  have  the  greatest  probability 
of  success  in  a  reasonably  short  time. 

3.1  Tunnel-Diode  Network 

Because  of  its  size  and  simplicity,  the  tunnel  diode  appeared  to  be  a  likely 
choice  for  an  active  element.  Properly  biased,  the  tunnel  diode  may  be  modelled 
by  a  negative  resistance  in  parallel  with  a  capacitance,  as  shown  in  Fig.  3-2. 

The  network  in  which  it  is  embedded  was  designed  on  an  intuitive  basis;  essentially, 
the  inductance  L  paralleling  the  diode  gives  a  pair  of  poles  in  the  right-half  plane, 
the  RLC  combination  gives  a  mirror-image  pair  in  the  left-half  plane,  and  Lq 
moves  the  poles  on  to  the  jio-axis. 
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Straightforward  network  analysis  yields 


z  =  4) 


4.  2/,VL 

P  P  [2LoLC  r2c2 


L  +2L 
o 

2  2 
L  L  C 
o 


If  Lq  is  chosen  to  satisfy 


L„C  2R2CZ 


RC  JLC1 


then  (3. 1)  becomes 


Z  =  (~) 


p(p2^) 


C  'n2  ■  _L 

P  LC 


RG 


2  Y2 

/lc/ 


For  p  =  ju .  Z  is  a  pure  reactance: 


2  2 
Uj{w  *u  ) 

z<jU)  -  JXU  •  i<t>  — Hr 


<“b"u  1 


where 


y  <  y 

a  b 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


Obviously,  X(y)  is  positive  for  0  <  y  <  y  and  negative  for  w  ^  w  ^  ~  and 

a  a  b 

thus  has  negative  slope  in  some  region  around  y  ;  this  is  the  desired  type  of 
behavior. 
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The  circuit  of  Fig.  3-2  was  constructed  and  tested.  The  predicted  frequency 
response  was  not  observed;  the  frequency  response  being  essentially  that  of  a 
passive  network.  It  is  felt  that  this  was  due  to  second-order  effects,  most  impor¬ 
tantly'  lead  inductance,  which  were  not  included  in  the  model  of  the  tunnel  diode. 
Further  effort  did  not  appear  justified,  and  the  circuit  was  abandoned. 


3.2  Negative  Impedance  Converter  Development 

Another  likely  active  element  for  synthesis  is  the  Negative  Impedance  Converter 
(NIC).  Its  use  is  readily  suggested  by  the  negative-reactance  characteristic  desired 
in  the  present  application;  in  addition,  a  number  of  active-synthesis  techniques 
using  the  NIC  have  appeared  in  the  literature. 

The  NIC  is  a  two-port,  active  network  with  the  property  that  its  input  impedance 
is  the  negative  of  its  load  impedance.  Any  linear  two-port,  active  or  passive,  may¬ 
be  characterized  by  its  hybrid  parameters: 


I 

+  h  _ 

v 

11 

1 

12 

2 

I 

■5"  h 

V 

21 

1 

22 

2 

(3.6) 

(3.7) 


If  the  two-port  is  terminated  at  port  2  by  an  impedance  Z  ,  then  the  input  im- 
pedance  at  port  1  is  given  by' 


Z .  =  h, 

m  11 


h12h21 


w 


ff  hU  =  h22  ~  aw*  X  hi2  h2i  =  l*  11160 


(3.9) 
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FIG.  3-3:  IDEALIZED  NEGATIVE  IMPEDANCE  CONVERTER. 
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FIG.  3-4a:  NEGATIVE  IMPEDANCE  CONVERTER 


!l  re  rb  R3 


FIG.  3-4b:  AC  EQUIVALENT  CIRCUIT. 
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FIG.  3-5:  SIMPLIFIED  HIGH-FREQUENCY 
MODEL  OF  NIC. 
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II 

o 

V 

h12  • 

21  ”  1  ’ 

h22  *  2»C  • 

(3.17) 

The  only  departure  from  ideal  behavior  now  is  the  non-zero  value  of  this  can 
either  be  absorbed  in  the  load  or  compensated  for  by  a  capacitance  2p  C  in  parallel 
with  the  input. 

The  final  version  of  the  circuit  is  given  in  Fig.  3-6.  Resistors  and  R^ 
are  bias  resistors;  since  at  signal  frequencies  they  appear  in  shunt  across  the 
input  and  output  respectively,  their  effect  is  cancelled  out  by  the  NIC  action. 

When  this  circuit  was  constructed,  some  difficulties  of  stability  and  reliability 
of  operation  were  encountered,  but  these  were  eventually  overcome.  The  bias 
resistors  must  be  adjusted  to  prevent  ’’lock-up”  (both  transistors  cut  off),  but  the 
adjustment  is  not  critical.  Oscillation  is  possible  for  certain  combinations  of  load 
impedance  and  generator  impedance  (impedance  seen  looking  to  the  left  and  the 
input).  For  most  measurements,  it  is  necessary  to  pad  the  input  with  a  resistance 
in  order  to  prevent  oscillation. 

Figure  3-7  shows  the  input  impedance  of  the  NIC  over  the  frequency  range 
1-5  MHz.  The  load  impedance  was  a  270  Q  resistor;  a  680 Q  resistor  was 
placed  in  series  with  the  input  during  measurement  to  maintain  stability,  but  its 
resistance  has  been  subtracted  from  the  measured  values  so  that  Fig.  3-7  gives 
the  actual  input  impedance  at  the  NIC  terminals. 

The  fact  that  the  negative  resistance  seen  at  the  input  is  on  the  order  of-600  O 
rather  than  -270 fJ  is  of  no  consequence,  since  this  is  merely  a  scale  factor  which 
can  easily  be  accommodated  by  changing  the  load. 
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Frequency,  MHz 

FIG.  3-7:  INPUT  IMPEDANCE  OF  NIC;  ZL  =  27CH-J0.  (*'s  indicate 
measured  values;  the  lines  are  merely  to  fuide  the  eye. } 
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The  rather  large  reactive  component  of  the  input  impedance  is  a  rather  more 
serious  problem,  it  indicates  that  the  ratio  of  input  impedance  to  load  impedance 
is  complex  rather  than  negative-real.  However,  the  circuit  is  still  in  the  early 
stages  of  development,  and  considerable  improvement  can  be  expected  as  the 
development  continues.  Furthermore,  the  NIC  need  not  be  perfect  in  order  to  be 
useful;  this  point  is  discussed  in  more  detail  below. 

Measurements  have  thus  far  been  restricted  to  the  1-5  MHz  frequency  band. 
This  is  sufficiently  high  for  the  effect  of  collector  capacitance  to  be  significant  and 
low  enough  that  the  frequency  variation  of  a  may  be  neglected,  and  is  thus  an  ideal 
range  for  testing  and  adjusting  the  compensation  scheme. 

Since  the  NIC  is  essentially  a  lowpass  device,  the  problem  of  increasing  tie 
bandwidth  and  the  problem  of  increasing  the  maximum  operating  frequency  are 
essentially  the  same.  The  critical  assumption  in  the  analysis  above  is  that  or  -  1, 
which  implies  that  operation  must  be  at  frequencies  well  below  alpha  cut-off.  To  a 
certain  extent,  then,  improvement  in  NIC  performance  depends  on  improvement  in 
transistor  performance.  Given  the  present  rate  of  development  of  transistor 
technology,  there  are  good  grounds  for  optimism  in  this  regard. 

There  is  also  the  possibility  of  compensating  for  the  frequency  dependence  of 
a  by  the  addition  of  more  passive  elements  to  the  circuit.  This  has  not  been  tried 
so  far,  since  compensation  for  collector  capacitance  is  a  much  more  pressing 
problem  in  the  lower  VHF  band,  but  there  is  no  reason  to  believe  that  it  cannot  be 
done. 

Even  if  perfect  compensation  cannot  be  obtained,  the  device  may  still  be  useful 
if  some  measure  of  NIC  action  takes  place.  Imperfections  in  the  NIC  can  be  over¬ 
come  by  changes  in  the  load  impedance.  As  a  trivial  example,  suppose  that  an 
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■ 


f 


[- 


input  impedance  of  - 1  -r  j  0  is  desired  at  some  particular  frequency.  Assume 
further  that  a  load  impedance  of  1  -rjO  produces  an  input  impedance  of  - 1  -  j  1 
rather  than  the  desired  value.  Provided  the  NIC  is  operating  linearly,  a  load 
impedance 


_  -1+jO 
L  -1-jl 


0.5  -  j 0.5 


should  give  the  desired  input  impedance. 

More  generally,  the  process  of  obtaining  a  particular  input-impedance 
characteristic  using  an  imperfect  NIC  would  require  determining,  as  a  function  of 
frequency,  the  load  impedance  needed  to  give  the  desired  input  impedance,  and  then 
finding  a  passive  network  realizing  (or  approximately  realizing)  the  necessary  load 
impedance. 

3.3  Conclusions 

The  most  important  result  of  the  experimental  program  has  been  the  develop¬ 
ment  and  preliminary  evaluation  of  an  NIC  with  frequency  compensation  for 
operation  in  the  VHF  region.  Since  NIC  development,  as  reported  in  the  technical 
literature,  has  generally  been  at  frequencies  below  10C  kHz,  the  results  which  have 
been  obtained  to  date  are  a  significant  improvement.  The  immediate  objective  of 
continued  development  in  this  area  is  to  obtain  usable  NIC  action  at  50  MHz  or 
higher.  This  objective  appears  quite  attainable. 
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Relationships  between  the  real  and  imaginary  parts  of  an  impedance  are  studied.  Given 
the  real  and  imaginary  parts  in  a  finite  frequency  baud,  a  technique  is  developed  for  finding 
these  parts  outside  the  baud.  Conditions  under  which  the  real  port  remains  positive  for  all 
u3  are  given. 


Use  of  a  negative  impedance  converter  for  realization  of  a  particular  impedance  arising 
from  impedance  loading  studies  is  discussed.  A  method  of  compensating  for  collector 
capacitance  and  extending  operation  of  a  transistorized  NIC  into  the  VII F  region  is  given. 


